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ABSTRACT 

Under the basic assumption that the observed turbulent motions in 
molecular clouds are Alfvenic waves or turbulence, we emphasize that the 
Doppler broadening of molecular line profiles directly measures the velocity 
amplitudes of the waves instead of the Alfven velocity. Assuming an equipartition 
between the kinetic energy and the Alfvenic magnetic energy, we further propose 
the hypothesis that observed standard scaling laws in molecular clouds imply a 
roughly scale-independent fluctuating magnetic field, which might be understood 
as a result of strong wave- wave interactions and subsequent energy cascade. 
We predict that oc is a more basic and robust relation in that it may 
approximately hold in any regions where the spatial energy density distribution 
is primarily determined by wave-wave interactions, including gravitationally 
unbound regions. We also discuss the fact that a scale-independent a% appears 
to contradict existing 1-D and 2-D computer simulations of MHD turbulence in 
molecular clouds. 

Subject headings: hydromagnetics — ISMxlouds — ISM: kinematics and dynamics 
— ISM: magnetic fields — turbulence — waves 



1. INTRODUCTION 

Motions of material within interstellar molecular clouds have long been known to be 
dominantly turbulent except in small low-mass dense cores. A significant advance in the 
study of turbulence and cloud support was the finding of a set of scaling relations, i.e.. 
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correlations among cloud quantities such as the characteristic radius R, mean non-thermal 
velocity dispersion a.^ and the average density n^,^ = p/m, for various samples of molecular 
structures ranging from 0.2 pc to 200 pc in size (Larson 1981; Leung et al 1982; Myers 1983; 
Dame et al 1986; Scoville & Sanders 1987; Solomon et al 1987; Falgarone, Puget & Perault 
1992; Blitz 1993), 

(T^ DC i?", a = 0.5, (1) 
njjj oc i?"^,7 = 1, (2) 

and 

(7^ oc n^f,/? = 0.5. (3) 

Only two of the above relations are independent; the power- law indexes are related by 
a = dj. Since the above values for a, (3 and 7 are controversial, we will refer Equations 
(l)-(3) as the Standard Scaling Laws. Largely established on a cloud-cloud basis (i.e., based 
on samples of clouds), the physical significance of these scaling laws lies in the natural 
extrapolation that they are on average a manifestation of self-similar individual clouds 
satisfying these relations (Larson 1981), although arguably these correlations can be due 
to observational selection effects or even artifacts (Kegel 1989; Scalo 1990; Myers 1991; 
Vazquez-Semadeni et al 1996). In general, observational determination of these scaling 
relations is based on average properties of clouds and subject to numerous uncertainties. 
Therefore these scaling laws are necessarily rough and perhaps only of value dimensionally 
in a strict sense. On the other hand, the large number of observational studies that do find 
more or less the same scaling relations indicate that these scaling laws may indeed tell us 
something physical (Scalo 1987; McKee 1989), at least in a dimensional sense. 

Proposed explanations to the scaling laws bifurcate into hydrodynamic and 
magnetohydrodynamic (MHD) regimes and were critically reviewed by Scalo (1987). Most 
of them converge on that one of the two independent scaling relations is due to the trend 
for virial equilibrium. The hydrodynamic approach takes the scaling laws as an evidence 
for energy cascade in strong eddy turbulence over a large range of scales, similar to the 
Kolmogorov law for incompressible turbulence (Larson 1981), which gives dimensionally 
(7^ ~ E}/^. In the MHD regime, however, comparably rigorous theoretical approach has 
been hampered by the lack of theories capable of describing strong compressible MHD 
turbulence. 



2. EVIDENCE FOR A SCALE-INDEPENDENT FLUCTUATING 

MAGNETIC FIELD 
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Given a relatively strong interstellar magnetic field and a moderate ionization fraction 
of 10"'' to 10~^ maintained by penetrating cosmic rays and UV photons, the gas in 
molecular clouds is expected to behave much like plasma which is subject to numerous 
instabilities. If one recognizes the existence of a general interstellar magnetic field and 
turbulence, then it seems difficult to reject the likelihood that the turbulent motions (even 
if it is initially purely hydrodynamic) further generate MHD waves such as Alfven waves. 
In fact, Alfven waves have long been proposed as a viable physical process for explaining 
the observed "turbulence" in molecular clouds (Arons and Max 1975, hereafter AM75; 
Zweibel & Josafatsson 1983). Identifying (implicitly or explicitly) the turbulent velocity as 
the Alfven velocity, it has been further proposed that the Standard Scaling Laws can be 
explained if kinetic energy, energy of the general magnetic field and gravitational energy are 
all in approximate equipartition and that magnetic field field B is largely a constant over 
the scales concerned (Pellegatti Franco et al 1985; Scalo 1987; Myers & Goodman 1988, 
hereafter MG88; Fleck 1988). The difficulty of understanding such a constant B, however, 
has been noted (Scalo 1987; MG88; Mouschovias & Psaltis 1995); it contradicts not only 
the "frozen-in" picture for the interstellar magnetic field (cf. Mouschovias 1987) but also 
the observational evidence ioi a. B — n scaling relationship (Troland & Heiles 1986; Heiles 
et al 1993). 

The Doppler broadening of the spectral line profiles, nevertheless, reflects the velocity 
amplitudes of the waves (i.e. bulk moving velocity of the particles), not the phase velocity 
Va- The velocity amplitude of an Alfven wave packet of scale ^ R is determined by the 
assumed equipartition between the kinetic energy density and the fluctuating magnetic 
energy density over the scale R (AM75) 



where 5v = {8ln2Y^^a^ and 5B = (8/n2)i/VB.[] Alfven (1953) was perhaps the flrst to apply 
this relation astrophysically (solar granulation). McKee & Zweibel (1995) and Zweibel 
& McKee (1995) have recently demonstrated this relationship for weak MHD turbulence 
(when 6B is much smaller than B). Although it is difficult to prove the validity of this 
relation in the strong compressible turbulence regime, both observational and theoretical 
studies of solar wind turbulence have confirmed this basic relation, which has become a 
basic building block of a lot of analyses (cf. Tu & Marsch 1995; Lau & Siregar 1996). Recent 



Gaussian distribution is assumed for the fluctuating magnetic field and the velocity field with root- 
mean-squares as and <Ty. Note also that all the physical quantities in Equation (4) are implicitly associated 
with scale R. 




— or 6v 



6B 



(4) 



(47rp)i/2' 
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MHD turbulence simulations for molecular clouds also verified this energy equipartition 
(Passot et al 1995; Gammie & Ostriker 1996). 

Next, let us look at the commonly assumed virial equilibrium for gravitationally bound 
clouds. Observationally, virial equilibrium is established by comparing the kinetic and 
the gravitational energy of gas (cf. Larson 1981; Bhtz 1993). Theoretically, however, the 

situation is much less trivial, as the presence of the general magnetic field should also be 
taken into account (cf. Nakano 1984; Mouschovias 1987; McKee et al 1995). Naively, the 
observationally established virial equilibrium seems to underscore the direct role of the 
general magnetic field, which is consistent with a possibility first raised by Mestel (1965) 
that the volume and surface terms due to the general magnetic field have opposite signs 
and thus tend to cancel. A potentially more interesting (and radical) possibility is that 
turbulent motions can effectively reduce the effective pressure or stress of the magnetic field 
(Kleeorin et al 1996). The gravitational energy density of a gas parcel at radius R from 
the center of gravity can be written as Uq = 'iCG ^'^^^ , where C is a numerical factor 
depending mainly on the density distribution and M{R) is the mass inside a radius it!. The 
simplest form for the virial equilibrium in accord with observations can be written as 

where m is mean mass of the molecules. This is the same approach as adopted by MG88 
and Caselli & Myers (1995). Defining a mean density pm so that M{R) — ^npmR^, the left 
hand side of Equation (5) becomes CGpmpR? ■ If Prn ~ P within the errors of observational 
measurements. Equations (4) & (5) lead to 

p = {-t:CG)-^'\^ + ^kTf/^R-K (6) 
3 47r m 



Equations (4) and (6) give 



It is clear that Equations (7), (6) and (4) match the Standard Scaling Laws, Equations 
(1), (2) and (3), respectively, if and only if o"^ is scale-independent and thermal energy 
is small compared to the non-thermal energy. Since the turbulent pressure in this case 
is Pturb — po'y = a scale-independent a% is consistent with the consensus that the 
scaling laws imply an approximate dynamic pressure equilibrium in molecular clouds 
(cf. Chieze 1987; Fleck 1988; Maloney 1988; Elmegreen 1989; Blitz 1993). For a sample 
of galactic molecular clouds associated with HII regions, Solomon et al (1987) found a 
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velocity dispersion-size relationship a^{km s~^) — 0.7(^)°-^. Prom Equation (7), we derive 
ub — jJL gauss — 24: fi gauss for C ~ 0.2 (uniform density). If these clouds are 

centrally condensed (C is somewhat larger), the required as will be correspondingly smaller. 
Note, however, when thermal kinetic energy density is not negligible on small scales, both 
the density-size and the velocity-size relations will be different from the standard forms. 
When thermal kinetic energy density dominates. Equation (6) yields the classic isothermal 
density distribution; otherwise the power-law index 7 for the density-size relation ought 
to be between 1 and 2, as discussed by MG88, McKee (1989), Fuller & Myers (1992) and 
CaseUi & Myers (1995). 



3. WAVE- WAVE INTERACTIONS: AN ANALOGY TO KOLMOGOROV 

CASCADE 

Since observations indicate that the turbulent kinetic energy is often comparable to 
the energy of the general background magnetic field in some molecular clouds (e.g., MG88; 
Crutcher et al 1993; 1994), MHD waves at different wavelengths are expected to have strong 
mutual interactions and thus become strong MHD turbulence (cf. Myers & Khersonsky 

1995) . Intuitively, non-linear wave- wave interactions and consequent energy cascade largely 
dominate the spatial energy density distribution (cf. Sagdeev & Galeev 1969; Elmegreen 
1990; Biskamp 1994; Goldreich & Sridhar 1995), similar to the Kolmogorov energy cascade 
for hydrodynamic eddy turbulence. Unfortunately, an eigenvalue description that is used for 
weak MHD turbulence breaks down in this case, leaving essentially no other effective means 
so far for establishing a working theory which could tell us what the spatial energy density 
distribution should be for strong compressible MHD turbulence. The theoretical difficulty 
is clearly witnessed by the recent controversy over the Iroshnikov-Kraichnan theory for 
incompressible MHD turbulence (e.g., Goldreich & Sridhar 1995; Ng & Bhattacharjee 

1996) , which is a much more simplified problem than strong compressible MHD turbulence 
in molecular clouds. 

In this case it seems useful to see what observations tell us. Specifically, as argued in the 
previous section, the observed scaling laws in molecular clouds indicate a scale-independent 
fluctuating magnetic energy density a%. In terms of energy spectrum, P^, i.e., energy per 
unit volume per unit wavenumber as a function of wavenumber /c ~ 1/i?, a scale-independent 
a% implies a power-law Pk oc k~^. What is intriguing is that this k"^ spectrum corresponds 
to equipartition of mode amplitudes (Biskamp 1994) or Raleigh- Jeans energy equipartition 
(Sagdeev & Galeev 1969) in the sense that the wave packets of different scales obtain an 
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equal share of the available energy within a unit volume of space. This implies further that 
wave packets on larger scales would dominate the total amount of energy. 

Now, is such a spatial energy density distribution physically plausible ? We consider 
two aspects of the problem. First, the energy flux of an Alfven wave packet on a scale 
R can be written as 5" = (cf. Lau & Siregar 1996). It is evident that an ensemble 

of interacting wave packets traveling in all directions would yield a net energy transport 
flux AS only in the direction of the negative spatial gradient for the magnetic pressure or 
energy density a^/An, and thus the developing tendency is to reduce the spatial gradients 
for 0"^. Second, a main characteristic of MHD turbulence is the existence of inverse cascade, 
i.e., energy cascade from small to large scales with the tendency to form self-organized 
large-scale coherent structures (Biskamp 1994), in contrast to the direct energy cascade 
(from large to small scales) in pure hydrodynamic turbulence. Therefore it is not unlikely 
that the interacting wave packets of different scales achieve a net tendency for an equal 
share of the available energy density due to the effective energy cascade in both the direct 
and inverse directions while the system is fully turbulent (cf. Biskamp 1994). 

One point that appears interesting to us is that if the kinetic energy does remain 
in equipartition with the fluctuating magnetic energy, then perhaps a hydrodynamic 
description and a MHD description could lead to comparable physical insights to the 
turbulence in molecular clouds, at least from an energy point of view. The MHD approach, 
however, does offer certain advantages. One of these is that MHD waves may effectively 
transport energy with a high speed and perhaps also a low dissipation rate without 
necessary bulky flow of gas, and yet achieve more-or-less a uniform pressure environment. 
The disadvantage of the MHD approach is, however, its mathematical complexity. 

Given that a scale-independent fluctuating magnetic fleld is physically plausible and 
the simplicity of it being able to explain the observed Standard Scaling Laws as discussed in 
previous section is remarkable, we feel obligated to propose this possibility as an empirical 
hypothesis for further observational and theoretical investigations. 



4. DISCUSSION 

4.1. (Tfi versus B 

A fluctuating magnetic field component is not only implied by the observed 
■turbulence" or hydromagnetic waves (Alfven 1953; AM75), but also expected to explain 
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why molecular clouds do not free-fall to the center of mass along the generally-ordered 
magnetic field lines (cf. Shu et al 1987). Shu (1991) even suggested the intriguing possibility 
that an anisotropic fluctuating magnetic field might explain the reported prolateness of 
low mass cores. The fiuctuating magnetic field as is, nevertheless, conceptually different 
from the general magnetic field B, even in the extreme case that cxb ~ B. This raises the 
issue as for if the Standard Scaling Laws require a constant general magnetic field B over 
a large range of scales, as some have discussed (e.g., Pellegatti Franco et al 1985; Scalo 
1987; MG88; Fleck 1988). The answer to this question is negative based on the above 
discussions, but a constant B would not contradict the scaling laws. A key point is that a 
constant B does not directly explain the Standard Scaling Laws without invoking as first, 
because the Alfven velocity is not the same as the observed fluctuating turbulent velocity 
amplitudes. In other words, even if S is a constant over a large range of scales, a separate 
physical mechanism would still have to be sought to produce a scale-independent as- In 
practice, however, it is likely that as becomes comparable to B on large scales. From the 
empirical B — n relation B ~ 1.5(^-^^73)^''^ gauss (Heiles et al 1993) and the estimated 




(7b ~ 24 gauss (see Section 2), it can be crudely estimated that as becomes comparable 
to B when — 260 cm~^, (T„ ~ 2.3 km s'^ and ~ 10 pc. Taking this at its face value, 
it seems that aB could be well below B, or equivalently (7„ considerably smaller than the 
Alfven velocity on scales smaller than a few parsecs. 



Generally speaking, our principal expectation is that regions where the spatial 
distribution of the turbulent energy density is dominantly determined by wave-wave 
interactions are more likely to demonstrate the Standard Scaling Laws, while regions where 
this is not the case are unlikely to do so. Perhaps this is why the recent observational 
results (Plume et al 1996) of massive star forming cores do not agree with the standard 
scaling laws, as localized energy sources might dominate the energy density distributions 
over the scales of massive star formation. Further, the velocity-density relation (Equation 
(3)) among the Standard Scaling Laws should be the most robust physical relation in the 
sense that it may be expected for any region where magnetic fluctuations become scale- 
independent. In contrast, the standard density-size and velocity-size relations require not 
only scale-independent magnetic field fiuctuations, but also a second physical relationship 




4.2. 



Predictions 
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such as the virial equihbrium^. This prediction seems supported by the observations of 
Falgarone, Puget & Perault (1992) which find a remarkable correlation cx ti^^^ for a 
sample of gravitationally unbound molecular structures, although it is possible that jump 
shocks may also give rise to this relation (J.Scalo, private communication). 



4.3. Comparison with Existing 1-D and 2-D MHD Simulations 



Our hypothesis for a scale-independent fluctuating magnetic fleld, however, appears 
to contradict the recent 1-D and 2-D computer simulations for strong compressible MHD 
turbulence in molecular clouds. These simulations report a fluctuating magnetic spectrum 
(energy per unit volume per unit wavenumber), oc /;;"*, where t ~ 2.0 — 2.3 for 1-D 
(Gammie & Ostriker 1996) and t ~ 1.85 for 2-D (Passot, Vazquez-Semadeni & Pouquet 
1995). This magnetic spectrum implies a fluctuating magnetic field, gb oc 
(John Scalo, private communication). In other words, these simulations show a steeper 
k-dependence for the magnetic spectrum than our hypothesis. Both simulations find 
essentially the same velocity spectrum (energy per unit mass per unit wavenumber) 
E'fc cx which is dimensionally consistent with the the observed standard velocity-size 
relation bv ~ B^l^ . It appears difficult for the simulations to explain the other observed 
Standard Scaling Laws. The difficulty of explaining the standard size-density relation 
in the 2-D simulations of Passot et al (1995), for example, is recently made clear by 
Vazquez-Semadeni et al (1996). These authors were then led to suggest that the standard 
velocity-size relation is the basic relation while the density-size relation may be an 
observational artifact, in clear contrast with our predictions in this paper. This raises the 
question as for if a simple dimensional analysis from a pure energy consideration as we have 
done in this paper is a valid approach for the MHD turbulence and if the equipartition 
between kinetic energy and fiuctuating magnetic energy is a good assumption. On the 
other hand, 1-D and 2-D MHD turbulence simulations in solar wind studies have been 
known to produce energy spectra too steep to compare with observations, possibly due to 
the restriction in dimensionality, enforced periodicity and other unrealistic simplifications 
(Marsch 1991; Tu & Marsch 1995). So it appears that the standard scaling laws and 
subsequently our hypothesis for a scale-independent fiuctuating magnetic field in molecular 
clouds remain a challenge to MHD simulations, or vice versa. We share the same hope with 



^Heithausen (1996) has recently reported a remarkable correlation n cx for high latitude clouds, 

which made him question if virial equilibrium is the true reason for a n cx R^^ relation. It is not clear yet, 
however, if his assuming equal distance for all the high latitude clouds would yield the correlations artificially. 
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the solar physicists that future 3-D MHD simulations will shed hght on this important issue 
in MHD turbulence. 

In conclusion, we proposed the hypothesis that strong wave-wave interactions have led 
to a scale-independent distribution of fluctuating magnetic energy density and consequently 
the observed scaling laws in molecular clouds. We suggest that further careful observational 
studies of interstellar atomic/molecular clouds not in virial equilibrium or not harboring 
active star forming activities may serve as valuable direct tests of the basic ideas presented 
in this article, especially in comparison to the contradicting predictions of recent 1-D and 
2-D computer simulations. 
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